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CN ' Abstract 

m . 

■ W bosons can be produced in the channels e^p + X ai HERA thus allowing 

to probe for anomalous trilinear couplings among the gauge bosons. We discuss the 
next-to-leading order (NLO) QCD corrections to the photoproduction of W bosons 
with finite transverse momentum at HERA. The higher-order QCD corrections re- 
, duce the factorization scale dependence significantly and modify the leading-order 

(LO) cross sections by ±0(10%). 
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■ 1 Introduction 

5t , The study of W bosons at different colliders serves as an important test of the Standard 

Model and possible extensions. W bosons can be produced at the ep collider HERA 
with a center of mass (cm.) energy ^/s ~ 318 GeV which is achieved by colliding elec- 
trons/positrons with energy Eg = 27.5 GeV and protons with energy Ep = 920 GeV. 
Since the production cross sections for the processes e^p — >■ e^W + X reach values of 
about 1 pb at HERA, the production mechanisms of W bosons can be studied and the 
existence of anomalous WlVj trilinear couplings can be probed [|lH^ . Moreover, W boson 
production represents an important SM background to several new physics searches. In 
particular it is the dominant SM process leading to isolated high energy lepton events 
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with missing transverse momentum |l|, ^]. In order to determine potential discrepancies 
between measurements and Standard Model (SM) predictions, the latter have to be suf- 
ficiently accurate and reliable. This is not guaranteed for the available LO calculations 
of W boson production H, ^. For an unambiguous test of anomalous contributions, 
it is necessary to extend the previous analyses to NLO accuracy. The first step in this 
direction has been made in Ref. 0, where the QCD corrections to the total resolved 
photoproduction cross section have been determined. However, the result cannot be used 
for W boson production with large transverse momentum which is dominated by direct 
photoproduction, the QCD corrections to which are discussed in this letter. 




= 7^ \ \ " log ^ - ^mix, I ^ - ^ ) J , (1) 




2 Leading Order 

W boson production at ep colliders is mediated by photon, Z and W exchange between 
the electron/positron and the hadronic currents of the process. In general two phase-space 
regions are distinguished: the deep inelastic (DIS) regime at large and the photopro- 
duction regime at small Q^, being the negative square of the transferred momentum 
from the electron/positron. The photoproduction cross section can be calculated by con- 
voluting the Weizsacker- Williams photon spectrum, 

\ ^ mm 

with the differential hadronic cross section with respect to the transverse momentum pxw 
and the rapidity yw of the W boson for 7g q'W: 

= 2pTw2^ -fl/e{x^)qp{Xp,fiF)— ^ T7 , 2 

apTwdyw ^ "'a;- S + Ui b atl 

where qp{xp, is the corresponding quark density of the proton with momentum fraction 
Xp and at the factorization scale fip- crtoill ~^ Wq') denotes the corresponding LO 
partonic cross section and x^ the photon momentum fraction of the electron/positron. 
(For the resolved part the Weizsacker- Williams spectrum /^/g has to be replaced by the 
convolution of f^/e with the corresponding parton densities in the photon and the partonic 
cross section by the resolved expressions.) The QED coupling at vanishing momentum 
transfer is denoted by a, the electron mass by rUe, and the minimal and maximal values 
of the photon virtuality by Qmin,Qmax- order to separate photoproduction from 
the DIS region we impose a cut of Q'^^.x ~ ^ GeV^ in the momentum transfer from 
the electron/positron line. The minimal value of is kinematically fixed, 

Qmin ~ '^e ^ _J ^ ) (3) 

where negligible higher order terms in the electron mass rrie have been omitted. The 
partonic Mandelstam variables are defined as 



ih + k 



2 



\2 



2 



t = {h-Pi)\ h = t-M^, 

u = (/c2-Pl)^ ui = u-M^, (4) 

where ki is the 4-momentum of the photon, k2 the incoming quark 4-momentum and pi 
the 4-momentum of the outgoing W boson. The momentum fraction Xp of the proton is 
given by Xp = x~ , where 

x,T, + 
~ x,S + U, ' 
and the lower bound of the x^ integration by 

X- = -^ll±^ (6) 

S,Ti,Ui denote the hadronic Mandelstam variables. They are related to the partonic 
variables as 

s = x.yXpS, ti = x^Ti, ui = XpUi. (7) 

The hadronic Mandelstam variables can be expressed in terms of the transverse momen- 
tum ptw and the rapidity yw (defined to be positive in the electron/positron direction) 
of the W boson, 

where |/o = | log Ep/ denotes the shift in rapidity between the laboratory system and 
the hadronic cm. system. 

The leading direct photon process •yq — > q'W (a typical contribution is depicted by the 
first diagram of Fig. 0) develops a singularity at LO when the final state quark q' becomes 
coUinear with the initial state photon. However, the finite transverse momentum p^w of 
the W boson has to be balanced by the final state quark, so that this singularity does not 
occur at LO for non- vanishing prw- The small region includes the contribution of the 
hadronic component of the photon giving rise to W+jet production via qq' — > Wg (second 
diagram of Fig. as a typical example) and the crossed processes gq{q) Wq'{q'). The 
treatment of the DIS region is straightforward (a typical contribution is shown in the 
third diagram of Fig. |l|). 

The result for the direct contribution can be cast into the form 

- Wq'] = _ GFM>ie,t,-ej,ur / , ^ ^ ^M^t) , (9) 
dti v2 s'^utf ^ ' 

where Gp denotes the Fermi constant, the W mass and e^/, tw the electric charges of 
the scattered quark g' and W boson, i.e. ew = eq — Cqi = ±1. The LO cross sections of the 
resolved processes qq' — ^ Wg and gq{q) — > Wq{q) are given by {Cp = 4/3 and Tr = 1/2) 



.[qq'-.Wg] = Z C 



dti ^ ^ 3a/2 sHu 
^[gqiq)^Wqiq)] = — , (10) 



3 



7 'VWWr 



Q Q 



K/WW W q 



fWTTW^ g e 
K/WW W q 



>7,Z 



K/WW W 



Figure 1: Typical diagrams of W boson production with finite transverse momentum at 
HERA: direct, resolved and DIS mechanism. 
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Figure 2: Dependence of the differential LO production cross sections for ptw = 40 
GeV on the cut Q'^^.x which separates the DIS and photoproduction regimes, after adding 
the DIS, direct and resolved contributions. The renormalization and factorization scales 
have been identified with the W mass, fj. = fj^p = fJ'R = M\y. 
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where denotes the strong couphng constant. 

The direct, resolved and DIS contributions add up to the total W ptw distribution. 
Direct photoproduction develops the dominant contribution, while the DIS part is smaller 
but competitive. The resolved component is negligible for ptw ^ 15 GeV 0. The 
consistency of the calculation requires that the dependence on the specific value of the 
cut Qmax which Separates the DIS and photoproduction regimes, should be small. This 
dependence is presented in Fig. ^ for the LO and W~ cross sections at pxw = 40 
GeV. The residual dependence is below the per-cent level and thus indeed sufficiently 
small. 



3 QCD Corrections 

For the dominant direct part we have evaluated the NLO QCD corrections. They consist 
of two parts, the virtual and real corrections. The virtual corrections are built up by all 
one-loop diagrams which are generated by virtual gluon exchange. Typical examples are 
displayed in Fig. ^ The real corrections originate from gluon radiation off the quark lines, 
see Fig. and the corresponding crossed contributions with the gluon in the initial state. 
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Figure 3: Typical diagrams of the virtual corrections to direct photoproduction ofW bosons 
at HERA. 
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Figure 4: Typical diagrams of the real corrections to direct photoproduction of W bosons 
at HERA. 



The Q CD-corrected result including virtual and real corrections can be cast into the 
form 



fl2 dir 



dpTwdy 
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with all parameters as defined in Eqs. (^-^ and hp denoting the corresponding parton 
density of the proton. The partonic virtual corrections are given by 

= ^S{s + h + m + M^) (12) 
dtidui dti 

which leads to Eq. with q replaced by amrt for LO kinematics. 
3.1 Virtual Corrections 

The virtual corrections have been computed via dimensional regularization in n = 4 — 2e 
dimensions. The quarks have been treated as massless particles. Due to the Ward- 
Takahashi identities no renormalization is required after adding all diagrams, since there 
is no input parameter at LO which is affected by QCD corrections. The required scalar 
loop integrals are taken from Ref. p[. The final double and single singularities of the 
virtual corrections arise from infrared gluon exchange and coUinear singularities due to 
the gluon and massless quarks, 

'-'virt^F 



dti TT dti 

L 

e2 e V2 ■ " ° -t 



C..rt = cJ-\-^-(l+log^]+Cy] (13) 



with 



and a lengthy finite term Cy. These singularities are cancelled by adding the real correc- 
tions and the counter terms due to the renormalization of the parton densities. 

As a cross check, the virtual corrections have also been evaluated by introducing small 
gluon and quark masses A, m. This computation has been worked out with the programs 
FeynArts FormCalc [1^ and Loop Tools [1^. The renormalization of the small quark 



mass has been performed in the MS scheme to obtain an ultraviolet-finite result. The 
renormalization, however, does not affect the result in the limit of the massless gluons 
and massless quarks, but has to be introduced to take the massless limit in a consistent 
way. The infrared and collinear singularities now appear as logarithms of the gluon and 
quark masses in the limit of small masses, 

CLrt = ^ ^ + ^ 3^ - 5 - ^ + ^4 

with a different finite term Cy. This result has to be translated into the MS scheme 
for zero masses. The corresponding piece to be added to the virtual corrections can be 



calculated from the results of Refs. |[TT|, |T^ and it is given by 



ACyirt = C, I log 3^ log — - 2 — + log + 2 ~ 
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e2 e V2 -t J 2 ^ -t 2 ^ -t j ^ ' 

After adding this transformation term to the virtual corrections C^j^ second ap- 

proach, both results are in agreement. 



3.2 Real Corrections 

The real corrections have been calculated by means of the massless dipole subtraction 



method introduced in Ref. |T^. The basic idea is to construct appropriate dipole terms 
^fjsub ^j^g^^ include all infrared and collinear singularities of the real matrix elements, but 
can be integrated out analytically up to the LO phase space. The final expression for the 
differential cross section is given by 

daNLO = d&Lo + [da'"'''^ - da'''''] + [da"'''' + d&l"^] + [daP"'' + c/^™''], (17) 

where da'^j^2 ^"^o parts of the dipole terms which were integrated out analytically. 

They cancel the infrared and collinear divergences of the virtual corrections (ia™^* and the 
collinear singularity of the counter term rfo"^"^* due to the renormalization of the parton 
densities at NLO, respectively. The NLO parton densities have been defined in the MS 
scheme. Each of the square brackets in eq. (p!7[) is individually finite. This procedure 
allows to calculate the real matrix elements in 4 dimensions. 

The dipole terms are constructed for each collinear and infrared configuration arising 
in the processes. They include the collinear singularities between final state partons 
themselves and those between the final state partons and the initial parton or photon. 
The collinear singularity related to the photon is absorbed by the renormalization of 
the direct part of the resolved photon densities. There is a subtlety in the numerical 
implementation of the dipole term corresponding to collinear 7 — splitting, see Fig. ^. 
Although this dipole term does not generate a new spurious singularity, the numerical 




Figure 5: Generic diagram developing the collinear singularity between the final state 
quark and the initial photon. 

integration turns out to be unstable, if the gluon in the final state is soft. In this case 
the quark emitted from the photon develops a large transverse momentum balancing the 
transverse momentum of the W boson. Since the dipole term has not been constructed 
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to cover this case numerically, we introduced an upper cut on the product of photon and 
quark momenta, 

2fciP2 < pit ■ (18) 

After introducing this cut into the dipole terms which are subtracted from the real ma- 
trix element, and the analytically integrated one the numerical integration is stable and 
independent of this cut. 



4 Results 

We analyze our final results for direct photoproduction of W bosons plus one jet W + jet 
as well as the inclusive process W + X, i.e. without defining jets. We use the inclusive 
kx algorithm for the jet definition with a cone size R < 1 |13|. Moreover, we require 



each detected parton in the final state to have transverse energy Et > 5 GeV. One-jet 
configurations arise in this framework from 3 different phase space configurations: (i) LO 
kinematics include only one parton in the final state and are thus one-jet contributions 
automatically; (ii) if both partons of the real corrections are inside the same cone of size 
R < 1, they are combined to one jet; (in) if one parton has Et < 5 GeV, it is undetected 
so that only the other parton is counted as a single jet. The contributions of these three 
configurations have to be added to obtain the cross section for VT-I- 1 jet production. The 
residual part of the real corrections is attributed to W + 2 jet production which has been 
added to the 1-jet part for the full W + X production rate. 



We have chosen CTEQ4M [14] and ACFGP [15| parton densities for the proton and 
the resolved photon, respectively. The strong coupling constant is taken at NLO with 
^MS ^ 202 MeV. At LO we used CTEQ4L parton densities with LO strong coupling 



and A5 = 181 MeV. In our numerical analysis we have chosen the following values for 
the input parameters: Mw = 80.43 GeV, Mz = 91.1876 GeV, a = 1/137.03599976, 
Gf = 1.16639 X 10"^ GeV-2, = 0.510998902 MeV. In the fermionic Z couplings 
we have introduced the Weinberg angle siia^ 9\y = 1 — M^/M"^, and the Cabibbo angle 
between the first two generations has been chosen as sin^c = 0.222 0. 

Setting = fip = Mw for the renormalization and factorization scales, we present 
the final results for the prw distributions of + 1 jet and W + X production in 
Fig. ^. The QCD corrections modify the direct contribution by about ±(10 — 15)% 
and they are thus of moderate size. To estimate the theoretical uncertainties, the renor- 
malization/factorization scale dependence of the direct contributions to the processes 
e^p + X is presented in Fig. ^ for HERA conditions. The scale dependence is 

significantly smaller, once the NLO corrections are included. The residual scale depen- 
dence is reduced from about 20% down to about 5%. Fig.|^ clearly indicates that the 
NLO QCD corrections are accidentially small at the central scale determined by the W 
boson mass. Since the uncertainties of the parton densities are of similar size, the total 
theoretical uncertainty can be estimated to be less than about 10%. 



^In our numerical analysis we neglected the contribution of initial b quarks. 
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Figure 7: Dependence of the direct contribution to W production on the renormalization 
and factorization scale = /x^ = /xr = ^Mw forprw — 40 GeV. The full curves represent 
the NLO predictions and the broken curves the LO scale dependences. The upper plot 
presents W + Ijet production and the lower the total sum ofW + X production. 
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In Fig. ^ we present our final results for the total transverse momentum distributions 
of and W~ production at HERA including the NLO QCD corrections to the direct 
photoproduction part. The individual contributions are also shown. The Figure clearly 
demonstrates the dominance of the direct photoproduction part, while the DIS part is 
smaller but still relevant. The resolved photoproduction part can safely be neglected for 
Ptw values larger than about 10 — 15 GeV. The rising total sum at ptw values below 
about 15 GeV signals the breakdown of the pure perturbative results and underlines 
the necessity of a soft gluon resummation to describe W boson production for smaller 
transverse momenta which however is beyond the scope of this work. 

The rapidity distribution of W boson production at HERA for ptw = 40 GeV is 
presented in Fig. | for all individual parts. The direct photoproduction part is shown at 
LO and NLO. As in the case of the transverse momentum distribution the QCD corrections 
to the direct contribution amount to less than about 10 — 15% and are of moderate size. 
They hardly change the shape of the rapidity distribution. It can clearly be inferred from 
this figure that there is a preference to produce W bosons at larger values of the rapidity 

yw- 

In order to allow a comparison of our results with existing analyses at HERA, we 
have determined the cross sections for transverse momenta below and above 25 GeV. 
In the first case we have used the NLO results of Refs. P, for the cross section with 
Ptw < 25 GeV, while for prw > 25 GeV we have integrated our results including the 
NLO corrections to direct photoproduction. We obtain the following values for the cross 
section in the case of a positron beam (the numbers in brackets are the LO values): 

Ptw < 25 GeV : a{W+) = 0.478 pb (0.363 pb) 

a{W-) = 0.484 pb (0.348 pb) 

Ptw > 25 GeV : a{W+) = 0.150 pb (0.143 pb) 

a{W-) = 0.084 pb (0.079 pb) . (19) 

No other cuts have been imposed. 



5 Conclusions 

We have presented predictions for W boson production at HERA including the QCD 
corrections to the dominant direct photon mechanism at finite transverse momentum of 
the W bosons. Working in the conventional MS scheme we find that the QCD correc- 
tions modify the direct contribution by about ±(10 — 15)% at the nominal renormaliza- 
tion/factorization scale fiR = fip = Mw and are thus of moderate size. In addition, the 
NLO corrections reduce the residual scale dependence of direct photoproduction to a level 
of about 5%. Taking into account also the uncertainties in the parton densities of the 
proton and the photon, the total theoretical uncertainty is estimated to be about 10%. 
However, the QCD corrections to the DIS part are still unknown. They are not expected 
to be significantly larger, because they have to cancel the Q'^ax dependence of the NLO 
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direct contribution and thus have to be of similar size. Since the QCD corrections are 
dominated by soft gluon effects, the shapes of the differential distributions are hardly 
affected. Therefore the results obtained in this work cannot explain the excess of isolated 
high-energy lepton events observed at the HI experiment Q]. 

Acknowledgements. 

We would like to thank S. Catani, C. Diaconu, R. Eichler, E. Elsen, M. Kuze, A. Mehta, 
P. Schleper and P.M. Zerwas for useful discussions. We are grateful to C. Diaconu, P. Schleper 
and P.M. Zerwas for carefully reading the manuscript. 



14 



References 

[1] J. Breitweg et al., ZEUS Collaboration, Phys. Lett. B471 (2000) 411. 

[2] U. Baur, J.A.M. Vermaseren and D. Zeppenfeld, Nucl. Phys. B375 (1992) 3. 

[3] M.N. Dubinin and H.S. Song, Phys. Rev. D57 (1998) 2927. 

[4] C. Adloff et al., HI Collaboration, Eur. Phys. J. C5 (1998) 575; 

HI Collaboration, Proceedings EPS2001, contributed paper 802, Budapest, 2001. 

[5] C.H. Llewellyn-Smith and B.H. Wiik, Report DESY 77/38; 
P. Salati and J.C. Wallet, Z. Phys. C16 (1982) 155; 
K. Neufeld, Z. Phys. C17 (1983) 145; 

A.N. Kamal, J.N. Ng and H.C. Lee , Phys. Rev. D24 (1984) 2842; 

G. AltareUi, G. Martinelh, B. Mele and R. Riickl, Nucl. Phys. B262 (1985) 204; 
E. Gabrielh, Mod. Phys. Lett. Al (1986) 465; 

M. B5hm and A. Rosado, Z. Phys. C34 (1987) 117 and Z. Phys. C39 (1988) 275; 
U. Baur and D. Zeppenfeld, Nucl. Phys. B325 (1989) 253; 

J. Bliimlein and G.A. Schuler, preprint PHE-90-21, in proceedings "Snowmass 1990, 

Research directions for the decade" ; 

M. Janssen, Z. Phys. C52 (1991) 165; 

C.S. Kim and W.J. Stirling, Z. Phys. C53 (1992) 601. 

[6] P. Nason, R. Riickl and M. Spira, Proceedings of the 3rd UK Phenomenology Work- 
shop on HERA Physics, Durham, 1998, J. Phys. G25 (1999) 1434. 

[7] M. Spira, Proceedings 'Workshop on Monte Carlo Generators for HERA Physics" , 
Hamburg, 1998, |hep-ph/9905469| . 

[8] K.Fabricius and I. Schmitt, Z. Phys. C3 (1979) 51; R.K.Ellis, D.A.Ross and 
A.E.Terrano, Nucl. Phys. B178 (1981) 421. 

[9] J. Kiiblbeck, M. Bohm and A. Denner, Comput. Phys. Commun. 60 (1990) 165; 

H. Eck and J. Kiiblbeck, Guide to FeynArts 1.0, Univ. Wiirzburg, 1992. 

[10] T. Hahn and M. Perez- Victoria, Comp. Phys. Comm. 118 (1999) 153. 

[11] B. Humpert and W.L.van Neerven, Nucl. Phys. B 194 (1981) 225; S.Dittmaier, 
Nucl. Phys. B565 (2000) 69. 

[12] S.Catani and M.H.Seymour, Nucl. Phys. B485 (1997) 291, (E) ibid B510 (1997) 
503. 

[13] S. Catani, Y.L. Dokshitzer, M.H. Seymour and B.R. Webber, Nucl. Phys. B406 (1993) 
187; 

S.D. Elhs and D.E. Soper, Phys. Rev. D48 (1993) 3160. 



15 



[14] H.L. Lai, J. Huston, S. Kuhlmann, F. Olness, J. Owens, D. Soper, W.K. Tung and 
H. Weerts, Phys. Rev. D55 (1997) 1280. 

[15] P. Aurenche, P. Chiappetta, M. Fontannaz, J.P. Guillet and E. Pilon, Z. Phys. C56 
(1992) 589. 



16 



